Elliptic flow at large transverse momenta from quark coalescence 
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We show that hadronization via quark coalescence enhances hadron elliptic flow at large p± 
relative to that of partons at the same transverse momentum. Therefore, compared to earlier results 
based on covariant parton transport theory, more moderate initial parton densities dN/dr](b = 0) ~ 
1500 — 3000 can explain the differential elliptic flow W2(px) data for An + Au reactions at ^/s = 130 
and 2QQA GeV from RHIC. In addition, V2(p±) could saturate at about 50% higher values for 
baryons than for mesons. If strange quarks have weaker flow than light quarks, hadron V2 at high 
p± decreases with relative strangeness content. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.Gz; 25.75.-q 



Introduction. The goal of relativistic heavy ion colli- 
sion experiments is to produce macroscopic amounts of 
deconfined partonic matter and study its collective be- 
havior. One of the important experimental probes of 
collective dynamics in A -1-^4 reactions is differential ellip- 
tic flow0, V2{pi.) = (cos(2(/)))p^ , the second Fourier mo- 
ment of the azimuthal momentum distribution for a given 
p_L . Measurements of elliptic flow at high transverse mo- 
mentum provide important constraints about the density 
and effective energy loss of partonsj^Q- 

Recent data from RHIC for Au-|-Au reactions at 
y^sjVTV = 130 and 200 GeV show a remarkable saturation 
property of elliptic flow in the region 2 GeV < p± < Q 
GeV with V2 reaching up to 0.2 [lUllll. The saturation 
pattern, which corresponds to a factor of two azimuthal 
angle asymmetry of high-p^ particle production relative 
to the reaction plane, is still waiting for theoretical ex- 
planation. 

The saturation and eventual decrease of V2 at high 
has been demonstrated as a consequence of finite inelastic 
parton energy lossQ. Though the qualitative features in 
the data were explained, for realistic diffuse nuclei the 
calculations show a rapid decrease of above p±^ > 3 — 4 
GeV contrary to the saturation out to « 6 GeV seen 
in the data. 

Calculations of elliptic flow based on ideal (nondissi- 
pative) hydrodynamics H, |^ 0, 0| can reproduce the 
low p_L < 2 GeV data at RHIC remarkably well, how- 
ever overshoot the data above > 2 GeV [j, [TJl ■ The 
lack of saturation is due to the assumption of zero mean 
free path and that local equilibrium can be maintained 
throughout the evolutionfslIT^. 

Covariant parton transport theorvj3. IT3 . IT5llT^IT7llT^ 
overcomes this problem via replacing the assumption of 
local equilibrium by that of a finite local mean free path 
\{s,x) = 1/ a{s)n{x). The theory then naturally interpo- 
lates between free streaming (A — 00) and ideal hydrody- 
namics (A = 0). Several studies confirm that initial par- 
ton densities d!J| and elastic 2—^2 parton cross sections 
estimated from perturbative QCD, dNg/dr]{h = 0) ~ 



agreement with the V2{pi_) data is possible, provided ini- 
tial parton densities and/or cross sections are enhanced 
by an order of magnitude to adNg/drj{h = 0) ~ 45000 
mb|3. A similar enhancement is indicated by the pion 
interferometry data as well^^- The origin of such an 
opaque parton environment is the RHIC "opacity puz- 
zle". 

To compare to the experiments, parton transport 
models also have to incorporate the hadronization pro- 
cess. The studies mentioned above considered two sim- 
ple schemes: Iparton Itt hadronization, motivated by 
parton-hadron duality, and independent fragmentation. 
An alternative model of hadronization is quark coales- 
cence, in which the relevant degrees of freedom are not 
free partons but massive (dressed) valence quarks. Glu- 
ons are assumed to have converted to quarks, therefore 
there are no dynamical gluons considered. 

Quark coalescence has been applied successfully in the 
ALCOrH^I and MICOrIIJ models to explain particle 
abundances and spectra in heavy-ion collisions. It was 
also suggested recently in Ref. |2J| as an explanation for 
the anomalous meson/baryon ratio and features of the 
elliptic flow data at RHIC. In this letter we show that 
hadronization via quark coalescence can resolve most of 
the "opacity puzzle" because it leads to an amplification 
of elliptic flow at high p± . 

Quark coalescence. The usual starting point of coales- 
cence models is the statement that the invariant spec- 
trum of produced particles is proportional to the product 
of the invariant spectra of constituents . This means 
that (assuming that different quark and anti-quark distri- 
butions are the same) the hadron spectra at midrapidity 
are given by those of partons via 



dNB 
d'^p± 

dNM 
d'p± 



iPi-) 



Cb{p±) 



CMiP±) 



dN. 



d^p 
' dN, 

_d'^p± 



(P±/2) 



(1) 



1000 and a-gg^ 
effects at RHIC 



3 mb, 



generate too small collective 
|. Nevertheless, quantitative 



where the coefficients Cm and Cb are the probabilities 
for qq meson and qqq baryon coalescence. We 
allow for pt dependent coalescence factors because more 
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careful treatment of the coalescence probleni[23j shows 
that such a dependence may arise, e.g., due to kmematic 
(energy) factors or strong radial flow. This, however, 
does not influence elliptic flow because it is a ratio from 
which the coalescence factors drop out (see Eq. Q). 

These relations are only valid for rare processes. This 
is not the case at high constituent phase space densi- 
ties, when most quarks recombine into hadrons and hence 
the number of hadrons is linearly proportional to that of 
quarks, dNh{p±) (x dNg{p±). 

At lower constituent densities, coalescence processes 
become relatively rare and therefore the usual coales- 
cence formalism works. On the other hand, most quarks 
hadronize via fragmentation into hadrons. Nevertheless, 
depending on how quickly the parton phase space den- 
sity drops with increasing p±, there can be a region of 
hadron transverse momenta that is populated dominantly 
via coalescence. The reason for this is that hadrons 
from coalescence have larger momenta than the aver- 
age quark momentum, dN^°°-\p±) = Ch[dNq{p±/n)]" , 
{n = 2,3), whereas hadrons from fragmentation carry 
only a fraction z < 1 oi the initial quark momentum, 
dNt'^'ip^) ^ dN,{pJz). 

At very low parton densities, e.g., at very high trans- 
verse momentum, the fragmentation process wins, in ac- 
cordance with the QCD factorization theorem. For exam- 
ple, a power law parton spectrum dNq/p±dp± ^ ApJ^ 
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Therefore, in heavy-ion collisions there can be three 
qualitatively different phase space regions. At very large 
transverse momenta particle productions is dominated by 
independent parton fragmentation. At lower transverse 
momenta coalescence prevails, which region can itself be 
subdivided into two parts: a very low p± (high phase 
space density) region where Eq. is not applicable, 
and a moderate density (higher p±) region, where Eq. 
(P) is valid. Because the density of produced particles 
depends on the centrality of the collision, the "bound- 
aries" of these regions de pen d on centrality. Only de- 
tailed quantitative studiesjSOl Isif of the relative contri- 
butions of the various hadronization processes, which is 
beyond the scope of this letter, could determine where 
the exact bounds are. Alternatively, the limits can be 
deduced from comparison with the experimental data. 

Anisotropic flow. For brevity we discuss only elliptic 
flow as the most important and interesting case. How- 
ever, Eqs. ©J Q, and all conclusions below also apply 
(i) when azimuthal anisotropies Vk{pj_) = (cos(fc(/)))pj^ of 
any order are present, and (ii) to any anisotropy coeffi- 
cient Vk instead of V2, even in the former most general 
case. 

In the coalescence region, meson and baryon elliptic 
flow arc given by that of partons via 

W2,Af(p±) ~ 2u2,g(^), V2.b{p±) ~ 3v2,q{^) , (2) 

as follows from Eq. and W2 ^ 1- For example, if par- 



tons have only elliptical anisotropy, i.e., dNq/p±dp±d^ 
{l/2TT)dNq/p±dp±[l + 2w2,gC0s(2$)], then 
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Fig. n illustrates the effect of quark coalescence on 
baryon and meson elliptic flow compared to parton ellip- 
tic flow. The latter is shown schematically by the solid 
line. At small transverse momenta, parton W2(p±) (x 
p\, as follows from general analyticity considerations. 
This region, before V2 becomes approximately linear in 
Pj^ could be relatively small (depending on the effec- 
tive mass of partons). At higher transverse momenta 
p_L > 1 — 2 GeV, parton elliptic ffow saturates as pre- 
dicted by parton transport Q, and then, possibly already 
above pj_ > 4 GeV, decreases according to pQCD parton 
energy loss calculations 0. The curve for baryon(meson) 
elliptic flow has been obtained by simply rescaling the 
parton curve by a factor three(two) both vertically and 
horizontally. As discussed above, for very low and very 
high p±, we boldly use Eq. 10) beyond its region of ap- 
plicability but doing so does not affect the discussion. 
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FIG. 1: Qualitative behavior of baryon and meson elliptic 
flow as a function of px from quark coalescence. 

There are three qualitatively different regimes in Fig.^ 

(i) In the small pj_ region where V2 (j>± ) increases faster 
than linearly, V2.b < V2,m < V2^q- It is not clear to 
what extent the coalescence picture is applicable in this 
region but it is interesting that the data does exhibit 
such a behavior. This ordering follows naturally from hy- 
drodynamics, where flow decreases with increasing par- 
ticle mass 9.,. . Id .llj . Similar mass dependence could also 
arise in a coalescence model because heavier hadrons can 
be formed by quarks with larger relative momentum (ig- 
nored in the current approach). 

(ii) In the intermediate pj_ region where W2(p±) depends 
linearly on transverse momentum, V2,b ~ V2.m- 

(iii) At large px, where parton V2{jpi_) increases slower 
than linearly, baryon flow becomes larger than that of 
mesons, V2.b > V2.m > V2,q, by as much as 50%. 
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Parton collective flow saturation, predicted for p± > 

1 — 2 GeV by parton transport 'S] j results in saturating 
meson/baryon flow at pj^ > 2 — 4 GeV that is amplified 
two/three-fold compared to that of partons. Saturation 
sets in at 50% higher p± for baryons than for mesons. In 
addition, any eventual decrease of parton elliptic flow at 
very high p±^ , would happen at two to three times larger 
p± for hadrons. 

The high-pj^ results above strongly differ from those 
obtained in Ref. ^2^. The reason is that, unlike Eq. (|2J), 
in Ref. 25j the coalescence of quarks was considered to 
be independent of their relative momenta and therefore 
hadron elliptic flow at high p±^ was similar to that of a 
high-px quark. 

If not all quarks show the same elliptic flow, further 
differentiation occurs because in that case 

V2,B=abc{p±) ~ W2,a(p_L/3) + W2,6(p±/3) + U2,c(p±/3) 
V2,M=ah{pi^) ~ W2,a(p±/2) + W2,6(p±/2) . (4) 

For example, strange quarks may have a smaller V2{p±) 
than light quarks, at high p± because heavy quarks are 
expected to lose less energy in nuclear medium j2y|, while 
at low p± due to the mass dependence of hydrodynamic 
flow. If v| < V2, elliptic flow decreases with increasing 
relative strangeness content within the baryon and meson 
bands, i.e., V2 > V2 ^ vf > V2 > v2 and > > ■ 
Possible solution to the opacity puzzle. While 
hadronization via Iparton In or independent frag- 
mentation approximately preserves elliptic flow at high 

2 < p_L < 6 GeV Q , quark coalescence enhances V2 two 
times for mesons and three times for baryons. Hence, 
the same hadron elliptic flow can be reached from 2 — 3 
times smaller parton V2, which requires smaller par- 
ton opacities, i.e., initial parton densities and/or cross 
sections. The amplification also allows the RHIC V2 
data to exceed geometric upper bounds derived based 
on a nuclear absorption model'2^ (the data are com- 
patible with those constraints only for idealistic "sharp 
sphere" nuclear distributions [2^ ) . Those bounds apply 
to the parton V2 and thus are two/three times higher for 
mesons /baryons. 

To determine the reduction of parton opacity quan- 
titatively, we rely on the results of Ref. [3 that com- 
puted gluon elliptic flow as a function of the transport 
opacity, x = / dzatrp{z) « adN / di] / {diO mb), from 
elastic parton transport theory for a minijet scenario 
of Au+Au at = 130.4 GeV at RHIC. Those re- 
sults can be conveniently parameterized as V2{pj_,x) — 
?;™'"^(x) tanh[pj_/po(x)]j where w™"^ is the saturation 
value of elliptic flow, while po is the p± scale above which 
saturation sets in. For the estimates here we assume that 
all gluons convert, e.g., via gg qq, to quarks of similar 
p± and hence ^Kpi) — (p^). 

As shown in Fig. |21 the increase of elliptic flow with 
opacity is weaker than linear, ~ x^'^^- There- 

fore, a 2 — 3 times smaller parton elliptic flow, which 
is needed to match the charged particle V2 data from 




X{b = 0) 

FIG. 2: Fit parameters 1;™°^ and po as a function of trans- 
port opacity x(b=0), for the V2{p±) results in Ref. 

RHIC in our coalescence scenario, corresponds to 3 — 
6 times smaller parton opacities adN/dri{b = 0) ~ 
7000 - 15000 mb than those found in Ref. 3]. The 
lower(upper) value applies if high-px hadrons are mostly 
baryons(mesons). Based on preliminary PHENIX data 
llil showing TTo/h^ « 0.5 between 2 < p^ < 9 GeV, 
one may expect mesons /baryons « 1, in which case 
adN/dri{b = 0) « 10000 mb. 

In Ref. only collective flow was considered and 
the parton opacity at RHIC was extracted using ellip- 
tic flow data from the reaction plane analysis. Tak- 
ing into account non-flow effects that contributed up 
to 15-20% I23 to the first elliptic flow measurements, 
parton opacities should be further reduced by 25% to 
adN/dri{b = 0) - 5000 - 10000 mb. For a typical elastic 
gg — > gg cross section of 3 mb, this corresponds to an 
initial parton density dN/d-q{b = 0) 1500 - 3000, only 
1.5 — 3 times above the EKRT perturbative estimate|19l|. 

The remaining much smaller discrepancy is compara- 
ble to theoretical uncertainties. For example, pertur- 
bative cross section and parton density estimates may 
be too low. If most hadrons formed via coalescence, 
the observed hadron multiplicity dNh/dr] « 1000 would 
imply much higher initial parton densities dN/dt] ^ 
2000—3000. Constituent quark cross sections, aqq « 4— 5 
mb, also point above the « 3 mb perturbative estimate. 
One effect that estimate ignores is the enhancement of 
parton cross sections a oc a^jp^ due to the decrease of 
the self-consistent Debye screening mass /i '-^ gTeffir) 
during the expansion. Finally, the contribution of inelas- 
tic processes, such as gg <-> ggg, to the opacity has also 
been neglected so far. A preliminary study shows'l6] 
that this contribution can be similar to that of clastic 
processes. 

Summary. In this letter we studied elliptic flow of 
hadrons formed from coalescence of quarks with simi- 
lar momenta. At high p± > 2 GeV we found an en- 
hancement of elliptic flow compared to that of partons. 
With the enhancement, moderate initial parton densities 
dNg/drj ^ 1500 — 3000 are sufficient to account for the 
charged particle elliptic fiow data from RHIC, providing 
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a possible solution to the RHIC "opacity puzzle" . At low 
p± < 1 GeV, on the other hand, hadron elliptic flow is 
suppressed. 

Quark coalescence gives a weaker baryon flow than 
meson flow at low p± < 0.5 — 1 GeV, while the op- 
posite, V2 > , at high p± > 2 — 3 GeV. Assum- 
ing all partons have similar elliptic flow, w l.Swl^ 
at high p±. If on the other hand strange quarks show 
weaker flow than light quarks, elliptic flow at high p± 
is ordered by relative strangeness content, such that 
> 4 « vf > > > vt, v^-"" > vf > i;f, 
and ^2 > f 2^ ~ /2. These predictions can be readily 
tested in current and future heavy-ion collision experi- 
ments. 

We emphasize that the quark coalescence picture and 
therefore our flow ordering predictions strongly rely on 
the assumption that quark degrees of freedom are domi- 
nant at hadronization. Therefore, experimental support 



for our predictions may indicate the formation of decon- 
fined nuclear matter in heavy ion collisions at RHIC en- 
ergies. 

We also note that at very high p± one expects a tran- 
sition from hadronization via quark coalescence to in- 
dependent fragmentation. An experimental signature of 
this may be the reduction of baryon V2 below meson V2 ■ 

When this work was in its final stage, two preprints 
addressing baryon to meson ratio at high p±, ||33| and 
[3ll |. were submitted to the arXiv.org e-print server. While 
these studies mainly focus on baryon and meson yields, 
the underlying physical arguments are very similar to 
those presented here. 
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